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Abstract 

The usual Cauchy matrix approach starts from a known plain wave factor vector r and 
known dressed Cauchy matrix M. In this paper we start from a matrix equation set with 
undetermined r and M. From the starting equation set we can build shift relations for some 
defined scalar functions and then derive lattice equations. The starting matrix equation set 
admits more choices for r and M and in the paper we give explicit formulae for all possible 
r and M. As applications, we get more solutions than usual multi-soliton solutions for 
many lattice equations including the lattice potential KdV equation, the lattice potential 
modified KdV equation, the lattice Schwarzian KdV equation, NQC equation and some 
lattice equations in ABS list. 
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1 Introduction 

In recent decades the research of integrability of difference equations has got remarkable progress. 
The property of multi-dimensional consistency [IH3] reveals integrability in some sense for dis- 
crete equations. Based on this property lattice equation defined on an elementary square can 
be classified [3] and the result is referred to as ABS list, which is surprisingly short and only 
consists of nine equations named as Q4, Q3, Q2, Ql, A2, Al, H3, H2 and HI. We list out these 
equations in Appendix IB"1 

ABS list has received a lot of attention. With regard to solutions, a lattice equation that is 
multi-dimensionally consistent provides automatically a Backlund transformation [3]. Then, by 
the Backlund transformations one can derive both seed solutions and 1-soliton solutions [4.5J. For 
V-soliton solutions, the relation between Q3 equation and NQC equation [6UT] can be used. NQC 
equation is derived in Cauchy matrix approach 0[8] and its solution can be expressed through 
Cauchy matrices. On the other hand, Q3 equation can be degenerated to the lower equations 
Q2, Ql, H3, H2 and HI. It then follows that A^-soliton solutions of these equations can be given 
in terms of Cauchy matrices [7]. Besides Cauchy matrix approach, several other approaches are 
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developed for finding multi-soliton solutions of the lattice equations in ABS list, such as bilinear 
approach [5], transformation and iteration method [S], Cauchy matrix approach with elliptic 
functions [10] . Inverse Scattering Transform [11U12] . and algebro-geometry approach [13^14] . 

The Cauchy matrix approach [7] (also see [5|I15][TH] ) . starts from a known plain wave fac- 
tor vector r and known dressed Cauchy matrix M, defines some scalar dependent variables, 
constructs their recurrence relations, and then from closed relations provides discrete lattice 
equations, such as the lattice potential KdV (lpKdV) equation, the lattice potential modified 
KdV (lpmKdV) equation, the lattice Schwarzian KdV (lSKdV) equation and NQC equation. 
Solutions of these obtained lattice equations can be given through the dressed Cauchy matrix. 

In this paper, instead of known plain wave factor vector r and known dressed Cauchy matrix 
M, we start from a matrix equation set consisting of three equations, among which two equations 
are used to determine plain wave factor vector r and the third equation is used to define M. We 
can build shift relations for some defined scalar functions and then derive lattice equations. This 
procedure we call generalized Cauchy matrix approach. In fact, the starting matrix equation set 
is demonstrated to admit more choices for r and M which leads to more solutions than usual 
solitons. In the paper we will give explicit forms of all possible solutions to the starting matrix 
equation set. As applications, these solutions are used to construct solutions for many lattice 
equations, such as the lpKdV, lpmKdV, lSKdV, NQC, Q3, Q2, Ql, H3, H2 and HI equation. 

The paper is organized as follows. In Sec. 2, we briefly review the Cauchy matrix approach. 
In Sec. 3 we describe the generalized Cauchy matrix approach. In Sec. 4 we solve the starting 
matrix equation set and in Sec. 5 as applications solutions for some lattice equations are given. 



2 Cauchy matrix approach 

As a preliminary part let us briefly review the Cauchy matrix approach. For more details one 
can refer to Ref. [7] or [16]. 

A Cauchy matrix is known as a square matrix G = (Gij)nxN, Gi,j = w l _ z . , where Zi ^ 
Wj £ C. We will use its symmetric form, which is 

G = (G id ) NxN , Gij = 7— j— r-, h + -kj € C. (2.1) 

The Cauchy matrix approach starts from a "dressed" Cauchy matrix 

M = (M iij ) NxN = FxGxH, M M = -^_ (2.2) 

where c$ G C, pi is called the plain wave factor defined as 

w= fE±*tW|±*L , )" rf (2.3) 



^p-hj \q-k 

with constants p, q, p®, and the dressing matrices are 

F = Diag(pi, • • • ,p N ), H = Diag(ci, • ■■ ,c N 

M satisfies the relation 



MK + KM = r *c, (2.4) 
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where 

K = Diag(fci, k 2 , ■ ■ ■ , k N ), r = (p 1 , p 2 , ■ ■ • , pn) T , *c = (ci, c 2 , • • • , c N ). (2.5) 

By~and^we respectively denote the shifts in n and m direction, i.e., f(n,m) = f(n + l,m), 
f(n, m) = f(n,m + 1). Then, from the basic shift relation 

p + ki _ q + h 

Pi = r- Pi, Pi = r^' 

p-ki q-ki 

one may build the following shift relations of M, 

M(pl + K) — {pi + K)M = r *c, (2.6a) 

(pi - K)M - M(pl -K)=r *c, (2.6b) 

M(ql + K) — (ql + K)M = r *c, (2.6c) 

(ql - K)M - M(ql - K) = r f C, (2.6d) 

where I is the N x N unit matrix. 

Next, for any (i,j) £ Z x Z, introduce a scalar function 

= * c Ki(i + M^KW, (2.7) 
which can be proved to have the symmetric property 

S {ij) = S UA. (2.8) 

It then follows from the dynamical relation (j2.6|) that one can reach to a set of recurrence 
relations: 

p S(ij) - 5W+1) = p 5(*J) + _ 5(*.o) l g'(Oj) j (2.9 a ) 

+ S (i,i+D = p 5(ij) _ g{i+ij) + S (i,o)g(o,j) t (2.9b) 

g 5(ij) _ = 9 5(*j) + 5(i+i j) _ 5(*,o) iS '(oj) } (2.9c) 

+ = g ^(ij) _ + 5(i,o)^(o,i)_ ( 2>gd ) 

Besides S^'^, there is another scalar function defined by 

5(a,&) = + K)' 1 ^ + M)- 1 (aI + K)~ 1 r, a,beC, (2.10) 
which also has symmetric property 

S(a,b) = S(b,a). (2.11) 

It can be proved that S(a,b) obeys the shift relations 

1 - (p + 6)5(a, 6) + (p - a)S(a, 6) = V{a)V{b), (2.12a) 

1 - (g + b)S(a, b) + (q- a)S(a, b) = V(a)V{b), (2.12b) 

where 

V{a) = 1 - fe(al + iTT^J + M") -1 r = 1 - *c(I + M)~ 1 (a/ + K) _1 r. (2.13) 

With the help of the crucial symmetric properties (j2.8j) and (|2. 1 1 1) . some lattice equations 
appear as closed forms of the recurrence relations (|2.9p and (|2.12|) . We list those equations in 
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Since S^' and S(a,b) are defined by the known elements M,K,r, *c, solutions of those 
lattice equations are therefore given apparently. 
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3 Generalized Cauchy matrix approach 



In the previous section we briefly introduced the Cauchy matrix approach. Following this ap- 
proach several lattice equations can be derived and their solutions are expressed through the 
known elements including the dressed Cauchy matrix M, plain wave factor vector r, constant 
diagonal matrix K and constant vector t c. 

In the following, we start with unknown M, r and K, and investigate a generalized Cauchy 
matrix approach. 



3.1 Recurrence relations 

Let us first give the following Lemma. 

Lemma 1. Suppose that matrices K,A 6 C^xN are anticommutative, i.e., 

KA + AK = 0, (3.1) 
where all the eigenvalues {k\,k2,--- ,k^} (some of them can be same) of K satisfy 

ki + kj^O, Vl<i,j<JV. (3.2) 

Then A is a zero matrix. 

We leave the proof in Appendix [Al Next let us start to derive some recurrence relations. 

Theorem 1. Consider the N x N matrices M = (Mij(n,m))N X N and K = (-fQj)7Vx7Vj and 
the N-th order vectors r = (pi(n, m), p2(n,m), ■ ■ ■ , pN(n, m)) T and l c = (ci, C2, • • ■ , cn), where 
Mij(n,m) and pj(n,m) are undetermined functions while Kij and Cj are constants. Suppose 
that M, r, K and *c obey the relations 

(pi - K)r = (pi + K)r, (3.3a) 
(qI-K)r = (qI + K)r, (3.3b) 
MK + KM = r f c, (3.3c) 

which we call the starting matrix equation set, where p, q are constants, I is the N x N unit 
matrix. Define scalar function 

S (i,j) = t c K j (I + M)~ 1 KV, i,j 6 Z, (3.4) 
then we have the following recurrence relations: 

p S^i) - <?W+i) = p s<M) + S {i+1 ' j) - S {i ' 0) S {0 ' j \ (3.5a) 

pS (ij) + S (W) = pS&i) - 5(i+lj) + S^S(°j\ (3.5b) 

q&ij) _ = q siiJ) + #(i+bi) _ S {i ' 0) S {0 ' j) , (3.5c) 

qS (i,j) + = q s(iJ) _ + gmigioj). (3. 5d ) 

Define another scalar function 

S(a,b) = t c(bI + Ky 1 (I + My 1 (aI + Ky 1 r, a,beC. (3.6) 
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If satisfies the symmetric property 

S&S) =sV> i \ (3.7) 

then S(a, b) is also of symmetric form 

S(a,b) = S(b,a) (3.8) 

and satisfies the shift relations: 

1 - (p + b)S(a, b) + {p- a)S(a, b) = V(a)V(b), (3.9a) 

1 - (q + b)S(a, b) + (q- a)S(a, b) = V(a)V(b), (3.9b) 

1 - (p + a)£(a, b) + (p - b)S(a, b) = V(b)V(a), (3.9c) 

1 - (g + a)S(a, b) + (q - b)S(a, b) = V(b)V(a), (3.9d) 

V(a) = 1 - *c(J + M)" 1 (al + K)~V. (3.10) 

ifere m i/ie following we suppose si ± K is invertible for s = 0,p, q, a, b, and all the eigenvalues 
{&i,&2,--- ,A;at} (some of them can be same) of K satisfy 

ki + kj^O, Vl<i,j<N. (3.11) 

Proof. Let us first derive a shift relation of K and M. From (|3.3cp and ()3.3ap and noting that 
(pi ± K)K = K(pl ± K), we have 

(pi - K)MK + K(pl - K)M = [pi + K)r *c. (3.12) 

Meanwhile, just left-multiplying pi + K on ()3.3cp yields 

(pi + K)MK + K(pl + K)M = {pi + K)r *c. (3.13) 

Subtracting (|3. 13|) from (|3.12p yields 

[{pi - K)M - (pi + K)M]K + K[(pl - K)M - (pi + K)M] = 0, 

which further, in the light of Lemma [TJ gives 

(pi - K)M = (pi + K)M . (3.14a) 

For q and hat-shift we have 

(ql — K)M = (ql + K)M. (3.14b) 

In the following, the proof is actually similar to the usual Cauchy matrix approach [7]. Taking 
~ shift of (|3.3cp we have 

r *c = MK + KM. 

and replacing the last term KM by ()3.14ap yields 

M(pl + K) - (pi + K)M = r *c. (3.15a) 
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Besides, if deleting the term KM from (j3.3c|) and (|3.14ap we have 

(pi — K)M — M{pl — K) = r *c. (3.15b) 

Similarly, using ()3.3cp and (|3.14bp we have 

M{ql + K) - {ql + K)M = f *c, (3.15c) 
(ql — K)M - M(ql — K) = r t c. (3.15d) 

Now, noting that the above shift relation ()3. 15[) is exactly the same as (|2.6p . and the matrices 
K, (al + K) and (bl + K)" 1 commute each other as being diagonals, one can get the recurrence 
relation (|3.5p as in [7]. For completeness, we show the procedure in the following. We introduce 
auxiliary vectors 

= (I + M)- 1 J RTV, (3.16) 
and S^'-?) defined in (13.41) is then rewritten as 



S (iJ) = t cK i(I + M)~ 1 KV = l cKiu^. (3.17) 
Using the shift relation (|3,3ap . from ()3.16p we have 

(I + M)fiW = K { r = {pi - K^K^pI + K)r, 

Le -> _ 

K\pl + K)r = {pi -K){I + M)u {i \ 

Then, employing the exchange relation 

{pi -K){I + M) = {I + M){pl -K) + r t c (3.18) 
that is indicated by (|3.15bp . one has 

K\pl + K)r = (1 + M){pl - K)u {i) + r 4 cm w , 
which further, left multiplied by (/ + M) , yields the relation 

{pi - K)u® = pu® + - 5 (i ' 0) ti (0) . (3.19a) 

Now replacing r in (|3.16p by using (|3.3ap . we have 

{I + M)u {i] = K r = {pi + KY 1 K i {pI - K)r, 

i.e., 

K\pl - K)r = {pi + K){I + M)u w . 
In this turn we make use of the exchange relation 

{pi + K){I + M) = {I + M){pl + K)-r *c 

that is derived from (|3.15ap . we will finally reach to 

{pi + K)u^ = puW - + 5 (i ' 0) S (0) . (3.19b) 



Symmetrically, we can derive shift relations for (q,~), 

(ql - K)u® = gtiW + - S^u^, (3.19c) 

(ql + K)u® = qu {i) - u {i+1) + 5 (i ' 0) M (0) . (3.19d) 

Now, with the shift relation (|3.19|) in hand, a left-multiplication of t cK 3 immediately yields the 
recurrence relation (|3.5[) for S^' 3 \ 

The symmetric property S^ 1 '^ = S^'^ plays an important role and we will discuss this 
property later. If this property holds, noting that the formal expansion 

it is easy to see that S(a,b) = 5(6, o), and similarly, for V(a) defined in (|3.10p . i.e., 

V(a) = 1 - l c(I + M)~ 1 (al + K)~ 1 r = 1 - l c(al + -K") _1 (J + M) _1 r, (3.21) 

the second equality of holds as well. 

In order to get the relation (|3.9p . we introduce an auxiliary vector 



u(a) = {I + M)- l (aI + K)- l r, (3.22) 

under which we have 

5(a, b) = *c(6I + K)~ 1 it(a), (3.23a) 
V(a) = 1 - *cu(a). (3.23b) 

Then, using (|3.3ap we have 

u(a) = (I + M) _1 (a/ + K)- X r = (I + M) _1 (al + K)- 1 (pI - Ky 1 (pl + K)r, 

Le -> _ 

(pi - + M)u(a) = r + (p - a)(al + K^r. 

Making use of (IBTTgj) and (|3.23bll yields 

(/ + M){pl - K)u(a) = (p- a)(al + K)~ 1 r + V{a)r. 

Then by left-multiplying (I + M)^ 1 we get 

{pi - K)u(a) = (p- a)u{a) + V»u (0) . (3.24) 

Next, left-multiplying *c(6J + K) _1 on (f!T24"P we can reach to 

1 - (p + 6)5(a, 6) + (p - a)S{a, b) = V(a)V(b), (3.25) 

i.e., the shift relation (|3.9ap . where we have used the expression (|3.23p and the symmetric form 
V(b) = 1 — *c(6J + K)~ l (I + M)~ 1 r which has been given in (j3.21|) . In a similar way we get 
relation (|3.9bp . ()3.9cp and ()3.9dp are derived from ()3.9ap and ()3.9bp thanks to the arbitrariness 
of a, b and the symmetric property 5(a, b) = S(b, a). 

Now the proof for Theorem [1] is completed. □ 
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3.2 List of equations 

With the assumption of the symmetric property (|3.7p . some lattice equations can be derived as 
closed forms of the recurrence relations given in Theorem Q] (see [7] for detailed derivation) . In 
the following we list them out together with their solution formulae. 

• Lattice potential KdV (lpKdV) equation: 

(p + q + w — w)(p — q + w — w) = p 2 — q 2 , (3.26a) 
w = 5(0.0) = t c y + M )-i r . (3.26b) 

• Lattice potential modified KdV (lpmKdV) equation: 

p{vv — vv) = q(vv — vv), (3.27a) 
v = 1 - 5(0-1) = 1 _ t cK -i(j + M yi r . (3 27b) 

• Lattice Schwarzian KdV (lSKdV) equation: 



(3.28a) 



[z — z)(z — T) q 2 

(z-z)(z-f) ~ P 2 ' 

z= t cK- 1 (I + My 1 K~ 1 r-z ---—, zq G C; (3.28b) 

p q 

• NQC equation (5]: 

1 - (p + b)g(q, b) + (p- a)S(a, b) _ 1 - (q + a)S(a, b) + (q - b)S(a, b) 

l-{q + b)S{a, b) + (q- a)S{a, b) 1 " (P + a ) 5 ( a > b ) + (P ~ b ) S ( a i h ) ' 

5(a, b) = l c{bl + K) _1 (J + M)~ l (al + K)~V, a, b G C. (3.29b) 

In the following our task is to solve the starting equation set (|3.3|) so that we can give explicit 
solutions of the above lattice equations. 



4 Solutions to the starting matrix equation set ( 13.31 ) 
4.1 Simplification and canonical forms 

In the starting matrix equation set (|3.3p l c = (ci, C2, • • • , cn) is a known constant vector and we 
need to look for the solution triad {K, r, M}. We suppose T is the canonical form of K, i.e., 

K = T"TT, (4.1) 

where T is the transform matrix. We define 

Mi = TMT- 1 , n = Tr, *c a = *cT _1 . (4.2) 
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It then follows from (|3.3p that 

(pi - r)ri = (pi + r)n, (4.3a) 
(gl-r)n = (gl + r)n, (4.3b) 

Mir + rMi = n fei. (4.3c) 

Besides, for the scalar functions S(a,b) and V(a) we find 

Si**) = t cK j (I + M)- 1 K i r= tir^I + M^^rVi, 
5(a,6) = *c(6/ + K)- 1 (/ + M)- 1 (a/ + K)~ 1 r = *ci(W + r) _1 (J + Mi)~ 1 (al + r) -1 ri, 
V(a) = 1 - *c(I + M) _1 (al + K) _1 r = 1 - *ci(I + Afi) -1 (al + r) -1 ri, 

i.e., -K" and its canonical form T lead to same S^', same S(a,b) and same V(a). One more 
invariant is 

Z(a,b) = < c(6I + K)- 2 (/ + M)- 1 (a/ + K)~ 1 r = *ci(fe/ + r)- 2 (/ + M 1 )- 1 (aI + r)' 1 ri, (4.4) 

which will be used in later discussion. Thanks to such invariance, in the following we neglect 
the subscripts in (j4.3j) and consider the simplified/canonical equation set 

(pi - r)r = (pi + T)r, (4.5a) 
(ql-Y)r = (qI + T)r, (4.5b) 
Mr + TM = r *c, (4.5c) 

where r is a N x TV matrix in canonical form. Here we would like to specify the invertible 
assumption for T. 

Proposition 1. Hereafter we always suppose T satisfies the following necessary invertible con- 
dition: si ±r is invertible for s = 0,p, q, a, b, and all the eigenvalues {k±, k2, • • ■ , &/v} (some of 
them can be same) ofT satisfy 

ki + kj^O, Vl<i,j<N. (4.6) 

We can then start from the elements in (|4.5p and replace Theorem Q] as the following. 

Theorem 2. Suppose that T is a N x N matrix in canonical form satisfying necessary invertible 
condition described in Proposition^ M, r, *c are defined as in Theorem{J\ and they satisfy the 
equation set (|4.5p . -By them we define 

S (i,j) = t cr i(/ + M)~ 1 rV, (4.7a) 

5(o,6) = l c(bl + r)" 1 (/ + MY x (aI + T)~ 1 r, (4.7b) 

V(a) = 1 - *c(I + M)" 1 (a/ + r) -1 r, (4.7c) 

Z(a, 6) = *c(6I + r)" 2 (/ + M)" 1 (al + T^r. (4.7d) 

Then, S^' obeys the recurrence relation (|3.5p . If S^' has the symmetric property 

s (i,j) = 5(j,0 j (4. 8 ) 

i/ieri S(a,b) and V(a) also have symmetric property 

S( a ,b) = S(b,a), (4.9a) 

V(a) = 1 - l c(I + MY l (aI + KY l r = 1 - *c(al + K)~ X (I + M)~ 1 r, (4.9b) 

and they satisfy the shift relation 
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4.2 Solutions to fl4~5]) 

Noting that in the equation set (j4.5|) *c is known and T is a matrix in canonical form, it is 
possible to give a complete discussion for the solution pair {r, M} of (|4.5p according to the 
eigenvalue structure of T. 



4.2.1 List of notations 

First we need to introduce some notations where usually the subscripts u and j correspond to 
the cases of T being diagonal and being of Jordan block, respectively. 



plain wave factor: pi 



p + h\ n ( q + h 



^p-hj \q-ki 
N-th order vector : r l ^({kj}i) = (pi,P2, ■■■ ,Pn) T 



p®, with constants p, q, p®, 



N-th order vector : r 



Pi 



N-th order vector 
iV-th order vector : PJ 

N-th order vector : g 



r[W] 



[N] 



' 1! 

(1,1, ,1)' 
(1,0,0,--- ,0) 5 



(iV-1)! 



a a z a 6 



N x N matrix : r^({%}f ) = Diag(Jfei, k 2 , 

(a ■■■ 
1 a ••■ 
1 a ■■■ 



N xN matrix : (a) 



N x N matrix 
N x N matrix 



\ • • • 

^(ft}f) =Diag( Pl , P2 , 
M JV1 ({c J }f) = Diag( Cl , C2 , 



• • , k N 

\ 




1 a J 
■■■ ,PN 
■■■ ,CN. 



( 



N x N matrix : Fj^ (fci) 



Pi 
2! 




Pi 
1! 






Pi 



\ 



iVxiV matrix: i3"^({ Ci }f) 



(iV-l)! 
/ Pi 

C2 
C3 

V CAT 



Pi 



(AT-2)! (7V-3)! 

■ ■ ■ Cat_ 2 CN-1 

■ ■ ■ CN-1 Cat 

c N 



\ 






P 1 J 

CN \ 




o J 



(4.10a) 
(4.10b) 

(4.10c) 

(4.10d) 
(4.10e) 

(4.10f) 
(4.10g) 

(4.10h) 

(4.10i) 
(4-lOj) 

(4.10k) 



(4.101) 
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N x N matrix : G l ^({kj}^) = (gij) NxN , g itj 



ki -\- kj 



Ni x N 2 matrix : G^j'^dkj}^ 1 ; a) = {gi,j)Ni*N 3 , 9i 



Ni. 



1 \3 



Ni x N2 matrix : G 



m' n * (o;6) = (gi,j)NixN 2 , gi,j = C]^^- 



h + a 



N x N matrix : G^a) = G^ ,JV] (a; a) = (gij)NxN, 9i,j = C 



a + &)*+■?- 1 ' 

1 (-iy +j 



+J'- 2 (2a) i +5- 1 ' 



(4.10m) 
(4.10n) 
(4.10o) 
(4.10p) 



where 



C' 



J 



3 i\(j-*y- 

The iVth-order matrix in the following form 



, (j>i)- 





( a 








• 





\ 




at 


«o 





• 







A = 


a-2 


ai 


a 


• 









\ ajv-i 


CLN-2 


aN-3 ■ 


• 01 


"0 


/ 



(4.11) 



NxN 



with scalar elements {cij} is a iVth-order lower triangular Toeplitz matrix. All such matrices 
compose a commutative set G w with respect to matrix multiplication and the subset 

G m ={A\Ae G lN] , \A\ ± 0} 

is an Abelian group. Such kind of matrices play useful roles in the expression of exact solution 
for soliton equations. For more properties of such matrices one can refer to Ref. 
For the notations we just listed out, it is easy to find the following facts. 

Proposition 2. 

(1) . Both Hj Nl \{cj}i) and G ! f\a) are symmetric matrices; 

(2) . Both T ! j\a) and Fj N \ki) belong to G !N1 and therefore they commute; 

(3) . For H ! r N1 {{ Cj }?) and VA € G [N1 , they satisfy 



m N \{ Cj }?)A=A T m N \{ Cj }?), 

or, in other words, Hj Nl ({cj}^)A is a symmetric matrix; 
(4)- When a / 0, T^a)" 1 is given by 



(4.12) 



rf(a)- 1 



1 

a 

-1 

— 








1 

a 

-1 




1 

a 


(_i)Vr-i 


(_l)Ar-a 


(_l)AT-3 




a"- 1 


a N ~ 2 



\ 





4 1 J 

a' a ' 



(4.13) 
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4.2.2 Basic cases 

Let us discuss solutions for the equation set (|4.5[) when T takes two basic forms. Notations can 
be referred to Sec 14. 2.1] 

Case i. r = r^({^}f). 

This case has been solved in Ref. [7J as well as in Secj2j Solutions to (|4.5p are 



r = r^{{k 3 }f), (4.14a) 

ki + kjJ Nxn' 



M = F^({k^)G^({k^)H^({ Cj }^ = (t^-).. ,.. (4.14b) 



In this case one gets soliton solutions [7J. 
Case 2. T = T [ J N] (k l ). 

This is also called Jordan block case. Motivated by the Jordan block case of many other 
soliton equations [18l[T9] , it is not difficult to find a solution for the equations (|4.5ap and (|4.5b|) , 

r = r l f ] (k 1 ). (4.15) 



To find a solution M to (|4.5cp . we first rewrite 

r [ f i (k 1 )=F-I l f i , f c = 1^-11, M = FGH, (4.16) 

where F = Fj^(ki), H = Hj^({cj}^) and G is an unknown matrix. It then follows from (|4.5c|) 
that 

FGHT + TFGH = F ■ if 1 • I l f T H, 
where T = T^(ki). Noting that the item (2) and (3) in Proposition [21 we then have 



which further reduces to 



FGT T H + FTGH = F ■ if 1 • J^ T H, 



GT T + TG = I\ N] ■I l f lT . (4.17) 



Thus we need to solve the above equation. To do that, we write 

G= (5i,52,- •• ,9n), (4.18) 
where {gj} are column vectors of G. Then from (|4.17|) we find 

(h + T)g j+1 + 9j =0, (j = 1, 2, • • • , N - 1). 

This indicates solutions 

g 1 =T [ J Nl (2k i y 1 ■I [ P=g^{2k x ), (4.19a) 

g j+ i = 9igm ^ la = 2kl , (j = l,2,-..,N-l), (4.19b) 

where we have made use of the item (4) in Proposition [21 Substituting (|4.19p into (|4.18p we 
can find 

G = G l f(h) 
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with explicit expression ( |4.10p[ ). 

Thus, in this case, the solution pair to (|4.5|) is given by 

r = r [ f(k 1 ), (4.20a) 
M = F^ih) ■ G^ih) ■ H™ ({ Ci }f ). (4.20b) 

Besides, using the commutative property of lower triangular Toeplitz matrices in the set G [7V1 we 
have the following result. 

Proposition 3. For an arbitrary constant matrix A € G lN! , Ar and AM solve the equation set 
(|4.5p , where r and M are defined in (|4,20D . 



4.2.3 Generic case 

To investigate solutions to the equation set (|4.5|) with a generic T, we need to consider some 
elementary equations. 

Lemma 2. Consider the following three elementary matrix equations: 

GT^Ukj}^) + T'^akj}^ )G = l' D Nl1 ■ I ! o llT , (4.21) 

GT l f 2l (bf + T^dkj} 1 ^ 1 ) G = I ! D Nl] ■ l'j N2lT , (4.22) 

G T ! j N2l (bf + T ! ? l] (a) G = If 17 • If 2]T , (4.23) 

where the unknown matrix G in the above three equations is of respectively Ni x N\, N\ x N 2 
and Ni x N 2 . For solutions we have G = G% 11 ({kj}" 1 ) solves (^2T|1 . G = G'of^dkj}^ 1 ; b) 
solves g22]), and G = G 1 ^ 2 ' \a;b) solves g23}. 

Proof G = G^dkj}? 1 ) satisfying (^2Tp can be verified directly. 
Let us consider equation (|4.22p . We suppose 

G = (gi,g 2 ,--- ,9N 2 ) (4.24) 

with column vectors {gj}- Then from ()4.22p we have 

(b + T^({k j }?))g 1 =I^, 

(b + T l ^ ] ({k j }?))g j+ i + g j = 0, (j = l,2,--. ,iV 2 -l), 
which further yields 

9j = (f^TTW feTW' • • • ' (kzkwT' C* = 1,2, • • • , n 2 ). 



The matrix (|4T24"P with above {gj} is just G [ £ ,Nl] ({fy}? 1 ; b) and it solves K22\i . 

The solving procedure for (|4.23|) is quite similar to the Case 2 in Sec 14. 2. 21 We rewrite G as 
(|4.24p and insert it into (|4.23p and then we find 



+ rj*> (a))g j+1 + 9j = 0, (j = 1, 2, • • • , N 2 
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which yields 



7\ ■ 



u 



These column vectors compose a matrix G = G l j N j 1:N2] (a;b) that solves <|£25|) . 



1,2,- •• ,iVo 



!)• 



Now we can investigate the case with a generic T: 

r = Dia g (r^ ll ({%}f 1 ),r! / iV2l (^ 1+ i),r! 7 JV3, (^ 1+2 ), 

For convenience we define column vectors 

/ r™({*,}*) \ 
rT 2] (k Nl+1 ) 
r l f 3] (k Nl+2 ) 



[N. 



(^JVi+(s-l)))- 



( ir \ 

rim 
1 j 

r[AT s ] 



Obviously, such r solves (|4.5aj) and (|4.5bj) with T defined in f|4.26[) and 



*c = (ci,c 2 , 
To solve (|4.5cp . we still suppose 



, CATj , Cat 1+ i, • • • , CjV 1 +jV 2 +-+A r s 



M = FGH, 

where 

F = Diag(Fr i ({^}f 1 ),^ JV2l (^ 1+ i),Fl JV3l (% 1+ 2), ■ ■ ■ .^(fcjVx+C-i))). 
^Diag^telf 1 ),^ 
G is a symmetric matrix with block structure 



G = G 



and each Gjj is a iVj x Nj matrix. Noting that 



(Gy) 



ft = i\f, at = r 7 if, 



from (|4.5cp we reach to 



GY T + TG = I G I T G . 



In terms of block structure it reads 

G 1 ^\{k^)+T^\{k i }f)G^ 



G ltj T^\k Nj _ 1+1 f + T^\{k j }^)G 1 , j = I™ ■ I [ ?' ]T , 



G i;j T l ^(k Nj _ 1+1 y + T^(k Ni _ 1+1 ) Gij = Pj 



, T 



[JV<] 



(K j < s), 
(K i < j < s). 
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Equations for Gj } \ and Gjj are the transpositions of (|4.31b|) and (|4.31cj) . respectively, due to 
G = G T . These equations in (|4.3ip are just the type of the three elementary matrix equations 
AMD, ([4T22]1 and P~23"]) . Thus based on Lemma [2] we have solutions 



Gi,i 
G 



G [ ^ ] ({k 3 }^;k Nj _ 1+1 ), (Kj<a), 



N 



h3 



G 



(k Ni _ 1+1 ,k N] _ 1+ i), (1 < i < j < s), 



and then from (|4.29p we get M. 

Let us conclude this part as the following. 



Theorem 3. For the equation set (|4.5p . i.e., 



(pi - T)r = (pi + T)r, 
(qI-T)r = (qI + T)r, 
MT + TM 



r *c, 



with generic 

r = Dia g (rr 1 ; 

and 



CATi+l),r^ 3y (A:Ar 1+ 2), 



'(fcjvi+(s-l))) 



(ci,C 2 , ••• , CjV!,CjVi+l, ■ •• , CjV 1 +jV 2 +-+ivJ, 



we /iaue solutions 



j 

Jn 3 1 



+1 



(kNi 
(^1+2) 



V rf'^fcjvi+^-i)) / 



M = FGH, 



where 



Diag« l7 ({^}f ), Ff»'(** 1+1 ), ^f 3 ^ 



,Fi^(fe JVl+(s _ 1) )) : 



H-Dia,g{H D {{Cjh ),Hj {{Cj} Nl+1 ),■■■, Hj {{Cj} Nl+N2+ ... +Ne _ 1+1 )), 
and G is a symmetric matrix with block structure 

G = G T = (G 

with 



G ltl = G^ ] ({k j }^), 
G ld = Gj A = G ! D N ^({k^;k Nj _ 1+1 ), 



Gij = Gj, = G l J N ; :N3l (kN t _ 1 +i,kN j _ 1 +i), 



(Kj<s), 
(Ki<j< s). 



Besides, for T, t c > r and M above mentioned, the pair 

{Ar, AM} 



(4.32a) 
(4.32b) 
(4.32c) 



(4.33a) 
(4.33b) 
(4.33c) 



(4.34) 
(4.35) 

(4.36a) 



(4.36b) 
(4.36c) 

(4.36d) 

(4.36e) 

(4.37) 
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also solve equation set (|4.5p with same T and l c, where 

A = Diag(I Nl ,A 2 ,A 3 ,- ■ ■ ,A S ), 

in which ijvi * s the N\-th order unit matrix and Aj is a Nj-th order constant lower triangular 
Toeplitz matrix. In fact, AT = TA. 

Here, {Ar, AM} gives a general solution to (|4.5p with a generic T. It is then quite easy to 
write out solutions for special cases. For example, when 

r= (TO(f.) r ™ ( ° Ki+i) ) . fe= («,«,-, (4.38) 
a solution pair to (|4.33p are 



r™(k Nl+1 ) J ' 



(4.39a) 



and 

;k Nl+1 



where 



M = FGH = F r ~, ^ T ,V \H, (4.39b) 

F= , F l o Nli ({^) \ fl _/flS u ({c,}f 1 ) 



Wru.Mfu o \ / JT™({c,}f; 

F™{k Nl+1 ) )> M W2] ({c,}S:f) 

A second example, when 



(4.39c) 



o iy 2] (fc 2 

a solution pair to (|4.33|) are 



and 



r=( J n V r[WaI ,, x ), fe = (ci,c 2 ,--- ,cjv 1+ iv 2 ), (4.40) 



(4.41a) 



r [ j m {k 2 ) 



G m {kl) G ^ N2 \ki-M) 
G^ N2] (h;k 2 f G^\k 2 ) 



M = FGH = F\ „ [n J 2] ,::'t ( 4 - 41b ) 



where 



F l /°\k 2 ) J' I Hf 2l ({c,}S+f) J" (441C) 



4.3 Symmetric property S'W) = 

The symmetric property S^'-*' = S^'^ directly leads to the symmetric property S(a, b) = S(b, a) 

and symmetric definition of V(a). Therefore this property is crucial for getting lattice equations 

from the recurrence relations (|3.5p and (|3.9p . In this subsection we will see such a property 
S (i,j) = S (j,i) holds for 

a generic T in canonical form. For the notations we use below please 

refer to Sec l4.2,ll 



16 



4.3.1 Two basic cases 

We start from two basic cases. 

Case 1. r = rf({k j }?). 

Solutions to (|4.5|) of this case are given in (|4.14|) . i.e., 

r = rf{{k j }?), M = FGH, (4.42) 

where 

F = FW({kj}?), H = Hf{{c^), G = Gf{{k j }?). (4.43) 

Also noting that 

r = F-lf, t c = I^ T -JT, 
from the definition (|4.7a|) we have 

sihi) = *cr J (i + M) _1 rv 

= lf T ■ HT j {I + FGH)~ 1 T i F ■ if. 

Since T, F and H are diagonal we can freely commute them and then have 

S (i,j) = z m T . ^((HF)- 1 + G)~ 1 r* • if. (4.44) 

Noting that is a scalar function and G = G T , taking transposition for (|4.44[) we immedi- 

ately reach to 

S (i,j) = s (i,j) T = r m T . + G)" 1 ^' • 1^ = S ij ' l \ (4.45) 

i.e., the symmetric property S^'i' = S^' l \ 

Case 2. T = Tf{k 1 ). 

In this case, the equation set (j4.5|) admits the following expression: 

r = rf{k l ) = F-lf, l c = lf T -H, M = FGH, (4.46a) 

where 

F = Ff{k l ), H = Hf({ Cj }?), G = Gf{k 1 ). (4.46b) 
Recalling Proposition [2] we know that 

G = G T , H = H T , FT = TF, HY = Y T H, HF = F T H, (4.47) 

in which the last equality indicates {HF) = (HF) T . With these in hand, inserting (|4.46p into 
the definition (|4.7ap we have 



S (i,j) = t c Yi (I + My^r 



lf T ■ HTi(I + FGH^rF ■ if 
W T . f-rTYiurtr . j?n vr\-l pri . rW {4.48) 



= if ■ (T T y H{I + FGH)" 1 FT 1 ■ I l J 
= lf T ■ (T T y {{HF)" 1 + G)" 1 ^ • If. 

Then, taking transposition we immediately reach to the symmetric property S^' = S^' l \ 
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4.3.2 Generic case 

Based on the above two basic cases, now for the generic case the discussion becomes easy. 
Go back to Sec. 14.2.31 For the equation set (|4.5p with generic 

r = Diag(r^ 1 H{%}f 1 ),r [ 7 iV2l (^ 1+ i),rr 3] (^ 1+2 ),-- - ,r™(fc JVl+(s _ 1) )) ) (4.49) 

one has 

r = F I G , t c = I G T H, M = FGH, (4.50) 

where I G is defined in (|4.27p . F, H and G are given in (|4,36p , A fact is that here T, F, H 
and G satisfy exactly the same relation (|4.47p . which means in this case will have a same 
expression as (|4.48p . and that means S^J') = S^' l \ 

Theorem 4. For the equation set (|4.5p with T in canonical form, £(w = t cT ] (I + M)~ 1 T i r 

has the symmetric property S^'^ = S^' 1 ' . 

5 Applications 

5.1 Solutions to some lattice equations 

Now, for the lattice equations we listed out in Sec l3.2l their solutions can be given as follows, 
respectively: for the lpKdV equation (|3.26a|) . 

w = S (0 ' 0) = *c(I + M)- l r; (5.1) 

for the lpmKdV equation (|3.27al) . 

t, = 1 - S^~ l) = 1 - fcT-^I + M) _1 r; (5.2) 
for the lSKdV equation (|3.28ap . 

z= t cT~ 1 {I + M)~ 1 T~ 1 r-z ---—, zq £ C; (5.3) 

P Q 

for the NQC equation (|3.29a|) 

S(a,b) = t c(6I + r)- 1 (I + M)- 1 (a/ + r)~ 1 r, a,6eC, (5.4) 
where r, *c, T and M are described as in Theorem [3l 

5.2 Solutions to the lattice equations in ABS List 
5.2.1 Re-parametrisation 

Ref. [7J builds the relation between NQC equation and Q3 equation in the ABS list as well as the 
relation between Q3 equation and the 'lower' equations Q2, Ql, H3, H2 and HI (also see |17j). 
In Ref. [7] these ABS lattice equations (see Appendix [Bj were re-parameterized so that their 
solutions can be expressed through Cauchy matrices. These re-parametrisations are [7J 

Q2,Q1: P=jpz^r, q = -qir^z, (55 \ 

uo . s. _ P _ 1 z _ Q _ 1 y ' J 

H2,Hi: p = -p 2 , q = -q 2 - 
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And the re-paramaterized lattice equations are 



Q3 : P(uu + uu) — Q{uu + uu) = (p 2 — q 2 ) I (uu + uu) + 



Q2 : (q 2 — a 2 )(u — u)(u — u) — (p 2 — a 2 )(u — u)(u — u) 
+ (P 2 -a 2 ){q 2 -a 2 ){q 2 -p 2 ){u + u + u + u) 
= (p 2 - a 2 )(q 2 - a 2 )(q 2 -p 2 )[(p 2 - a 2 ) 2 + (q 2 - a 2 ) 2 - (p 2 - a 2 ){q 2 - a 2 ) 



Ql : (q 2 — a 2 )(u — u)(u — u) — (p 2 — a 2 )(u — u)(u — u) 



8 2 a\p 2 -q 2 ) 



H2 
HI 



(u — u)(u — u) = p 2 



(p 2 — a 2 )(q 2 — a 2 ) ' 
H3 : P(a 2 - q 2 )(uu + uu) - Q(a 2 - p 2 )(uu + uu) = 5(p 2 - q 2 ), 
(u — u)(u — u) + (p 2 — q 2 )(u + u + u + u) = p 4 — g 4 , 

I 2 

where in (|5.6a|) (p, P) = p and (q, Q) = q are the points on the elliptic curve 

{(x,X)\X 2 = (x 2 -a 2 )(x 2 -b 2 )}, 

in (IfTfidl) 

P 2 = a 2 -p 2 , Q 2 = a 2 -q 2 , 
and in Q3 and Q2 the dependent variable u has been scaled by 

u — > u(b 2 — a 2 ), u 

respectively. 
5.2.2 Solutions 



4 

a u 



(p 2 — a 2 ) 2 (q 2 — a 2 ) 2 ' 



(5.6a) 



(5.6b) 

(5.6c) 

(5.6d) 
(5.6e) 
(5.6f) 

(5.7) 



In Ref. [7] solutions to Q3 equation (|5.6ap is given by means of the relation between NQC 
equation and Q3 equation; solutions for other equations in the list ()5.6[) are derived by using 
the degeneration relation 



Q3 



Q2 



Ql 



H3 




H2 




HI 







Fig.l Degeneration relation 

In the following we list out solution formulae [7] for the lattice equations listed in (|5.6p . and 
here we claim that these formulae are valid for all the solutions {r, M} of our starting matrix 
equation set (|4.5p together with corresponding T and 'c which are described in Theorem [3l 

For Q3 equation (|5.6aD . 

u = A • G(a, b) + B- G{a, -b) + C ■ G{-a, b) + D ■ G(-a, -b), (5.8a) 
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where 

G <* »> = ( (p - a)( P - b) ) " ( {q - a )[q - b ) ) "P - <° + 6 ' S <°' < 5 *> 
S(a,b) is defined in (|4.7bj) . and A,B,C,D are arbitrary constants satisfying 

AZ)(a + 6) 2 -BC(a-b) 2 = (5.8c) 

16ao 

For Q2 equation (|5.6bp . 

u =7(0) [-(£ 2 + 1) + a£S(-a, a) + a 2 (Z(a, -a) + Z(-a, a)) 

4 n -i 

+ AD + — p(a)(l - 205(0, a)) + — p(-a)(l + 2aS(-a, -a))J , (5.9a) 

where 

/ pn om \ 

4 = 2a ^ ^ + 7T + £ > (4o a constant), (5.9b) 

\a z — p z a z — q z / 

/p + a\ n / g + a\ m 

P( a ) = Poo, ( Poo a constant), (5.9c) 

\p-aj \q-aj 

a 4 

^ = (p2_ a 2)2 ((? 2_ a 2)2> ( 5 - 9d ) 

S(a,b) and Z(a,b) are defined in (|4.7p . and A, D are arbitrary constants. 
For Ql equation (|5.6c|) . solution with exponential background is 

u = Ap{a){\ - 2aS(a,a)) + D p{-a)(l + 2aS(-a, -a)) + B(£ + 2aS(-a,a)), (5.10a) 

where S(a,b) is defined in (|4.7bp . £ defined in (|5.9bp . p(a) defined in (|5.9cj) . and A,B,D are 
arbitrary constants satisfying 

16AD + AB 2 = 5 2 ; (5.10b) 
while solution with linear background is 

u = 8{v 2 - 2vS^~ 1) + 25 ( " 1 '- 2) ) + 2A(y + c - S^ 1 '^), (5.11a) 

where is defined in (|4.7ap . 

Tl Tfl 

v = 1 h Vq, (uq a constant), (5.11b) 

and A, cq are arbitrary constants. 
For H3 equation (|5.6d|) . solution is 

u = (A + {-l) n+m B)W{a) + (D + (-l) n+m C)^- 1 T/(-a), (5.12a) 

where V{a) is defined in (|4.7c|) . -d is defined as 

■d = (-^—) ( tfoo, ( #oo a constant), (5.12b) 
\a-pj \a-qj 
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and A,B,C,D are arbitrary constants satisfying 

r 

AD - BC = —. (5.12c) 
4a 

For H2 equation (j5.6e|) . solution is 

u = ic 2 - C5 {0 ' 0) + 2S(°' 1 ) - ^ 2 + (-l) n+ "M(C + c - 2S , (°'°)), (5.13a) 

where S^*'" is denned in ()4.7ap . £ is defined as 

C = pn + + Co, ( Co a constant), (5.13b) 

and A, cq are arbitrary constants. 
For HI equation ()5.6f|) . solution is 

u = A(( - SW) + {-l) n+m B{( + c - 2S(°'°)), (5.14a) 

where is defined in p~7aj) . C defined in (|5.13bp . A, 5, Co are arbitrary constants and A, B 
satisfy 

,4 2 - = 1. (5.14b) 
Here we note that by the transformation 

w = C — u, 

the lpKdV equation (|3.26ap is transformed to HI equation (|5.6f|) . This corresponds to A = 
1,B = in (|5.14p . Besides, we do not consider Al and A2 equation in the ABS list due to 
simple relations with Ql and Q3(<5 = 0) equation. 

6 Conclusions 

The usual Cauchy matrix approach starts from known plain wave factor vector r and known 
dressed Cauchy matrix M. In the paper we have described a generalized Cauchy matrix ap- 
proach where the starting point is a matrix equation set with unknown r and M. Such a 
starting equation set can equally lead to those recurrence relations for the defined scalar func- 
tions (SW) and S(a,b). With regard to solutions, we showed that K in the starting equation 
set provides same solutions for the related lattice equations as its canonical form T dose. This 
enables us to simplify the starting equation set and investigate its solutions according to the 
eigenvalue structure of the canonical matrix Y. We obtained explicit forms of r and M for a 
generic canonical matrix V. Besides, we proved the symmetric property of S^w = S^'*', which 
is crucial to generating closed recurrence relations (lattice equations). As applications, based on 
Ref. [7] we obtained solution formulae for many lattice equations, such as the lpKdV, lpmKdV, 
lSKdV, NQC, Q3, Q2, Ql, H3, H2 and HI equation. Since the starting equation set admits 
more choices for r and M, solutions of the above lattice equations are more general than usual 
solitons. In fact, in some sense, solutions corresponding to T being a Jordan block can be con- 
sidered as a limit result of the case of T being diagonal (cf . [T81 - I20] ) , or multiple-pole solutions 
in direct linearization approach. Besides, usually, zero eigenvalues lead to rational solutions |21j . 
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However, since we need T to satisfy invertible conditions (see Proposition [TJ) , eigenvalues of T 
can not be zero and therefore here the obtained solutions do not include rational solutions. 

Finally, let us go back to the starting equation set (|4.3p or (|4,5p . which consists of three 
equations. Reviewing the generalized Cauchy matrix approach, we can see that the first two 
equations are used to determine the plain wave factor vector r. Actually, these two equations 
include implicitly 2nd-order dispersion relation, which determines the obtained lattice equations 
are of the KdV-type (cf. |7]). If the dispersion relation is of 3rd-order, one can get lattice 
Boussinesq type equations j22H25j . The third equation in the starting equation set is used to 
define M. To find the explicit form of M we write it in the dressed form M = FGH. This 
works as well as in the proof of the symmetric property S™' = S^' l \ Besides, some algebraic 
properties and technique also make the discussions and proofs neat. 
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A Proof of Lemma [L] 

There is a nonsingular matrix T € Cn x n such that 

r = TKT^ 1 , (1.1) 

where T = ("fij)NxN is an upper triangular matrix with = ki for i = 1,2, •• • , N. We also 
denote 

B = TAT' 1 = (/3i, ,p N ), 
where {(3j} are column vectors of B. Then (|3.ip can be rewritten as 

TB + BT = 0. (1.2) 

The first column of (|1.2p reads 

(feif +r)/3i = o. 
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Noting that \k±I + T\ ^ due to the condition (|3.2|) . we have (3i = 0. Then, with this in hand 
the second column of (|1.2|) is simplified to 



which yields @2 = 0. Repeating this procedure step by step we can successively get (3% = (3^ = 
■ ■ ■ = /3tv = 0. This means B = and then A = 0. 

B ABS lattice equations 

Here we list out lattice equations in ABS list [3]. These equations are 
QA : p(uu + uu) — q{uu + uu) 



(k 2 i + r)/3 2 = o, 




7j ((uu + uu) — pq(l + uuuuj) , 



(2.1a) 



Q3 : 



p(l 



q 2 )(uu + uu) — q(l — p 2 )(uu + uu) 



V ^pq J 

p(u — u) (u — u) — q(u — u) (u — u) + pq[p — q) (u + u + u + u) 

= M(p-q){p 2 -pq + 4 2 ), 




(2.1b) 



Q2 : 



H2 



HI 



H3 



Ql 



Al 



A2 



p(u — u)(u — u) — q(u — u)(u — u) = S 2 pq(q — p), 
p(uu + uu) — q(uu + uu) = 5(q 2 — p 2 ), 
(u — u)(u — u) = (p — q)(u + u + u + u) + p 2 — q 2 , 
(u — u)(u — u) = p — q, 

p(l — q 2 )(uu + uu) — q(l — p 2 )(uu + uu) — (p 2 — q 2 )(l + uuuu) = 
p(u + u)(u + u) — q(u + u)(u + u) = S 2 pq(p — q), 



(2.1c) 
(2.1d) 
(2.1e) 
(2-lf) 
(2-lg) 
(2.1h) 
(2-li) 



where p, q are spacing parameters, 5 is an arbitrary constant, and in Q4 equation (p,P) and 
(g, Q) are points on the elliptic curve [T7] 

{(x,X)\X 2 = x 4 -jx 2 + 1}. 
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